function [var_new] = limitRate(var, varp, h, varRateLimit)
   if (var - varp)/h > varRateLimit
       var_new = varp + varRateLimit*h;
   elseif (var - varp)/h < -varRateLimit
       var_new = varp - varRateLimit*h;
   else
       var_new = var;
   end
end